Measurement of the rare decay 7r° —> e + e 
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The branching ratio of the rare decay tt° — » e + e~ has been measured precisely, using the 
complete data set from the KTeV E799-II experiment at Fermilab. We observe 794 candidate 
7T° — > e + e~ events using Kl — > 37r° as a source of tagged 7r°s. The expected background is 
52.7 ± 11.2 events, predominantly from high e + e~ mass ir° — > e + e~7 decays. We have measured 
B(tt° -> e+e", (m e+e -/m^o) 2 > 0.95) = (6.44 ± 0.25 stat ± 0.22 syst ) x 10" 8 , which is above the 
unitary bound from 7r° — > 77 and within the range of theoretical expectations from the standard 
model. 



I. INTRODUCTION 

In this report we present a new measurement of the 
7T° — > e + e~ branching ratio using a larger data set from 
KTeV-E799 at Fermilab. This result used all data taken 
in the two runs of the experiment (1997 and 1999-2000). 
It supersedes the previously published measurement |]J 
from KTeV-E799, which used only the 1997 data. The 
basic measurement technique of using Kl —> 37r° as a 
source of tagged ir° decays is adapted from the previous 
analysis. 



A. The Decay tt° -> e + e~ 

The rare decay tt° — > e + e~ proceeds, to lowest order, 
in a one-loop process via a two-photon intermediate state. 
The decay rate was first predicted by Drell and has 
since received considerable attention both theoretically 
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and experimentally. Relative to the ir° — > 77 rate, it is 
suppressed by two powers of a and is further suppressed 
by 2{m e /m 7r o) 2 due to the approximate helicity conserva- 
tion of the interaction. The lowest order contribution has 
been calculated exactly in terms of a form factor Q , and 
lowest order radiative corrections have been calculated 
The contribution to the rate from on-shell photons 
is model independent and can be calculated exactly to 
form a lower "unitary bound" [j| on the branching ratio, 
B(7T° — > e + e~)> 4.69 x 10 -8 , neglecting radiative correc- 
tions. 

The primary interest in the decay rate is the excess 
above the unitary bound, as this is the contribution from 
virtual photons. Attempts to model the form factor and 
make predictions for the off-shell photon contribution 
have been made, most successfully using vector meson 
dominance (VMD) and chiral perturbation theory (%PT) 
approaches [a, la, 13 • A new measurement is significant 
for xPT, where 7r° — > e + e~ represents a tight experimen- 
tal constraint on calculations. It is of particular interest 
because 7r° — > e + e _ is the best-measured decay of a pseu- 
doscalar meson to a lepton pair and has no significant 
contributions from short-distance physics. Constraints 
on xPT from 7r° — ► e + e~ can be used to improve predic- 
tions for other P° — > l + l~ decays, including r\ — > 
and the long-distance contribution to — > fi + yT . The 
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smaller short-distance contribution to K a L — ► /i + ^t~ is 
dominated by a top-quark loop and thus is a potential 
source of information on |Vtd| if the long-distance contri- 
bution can be subtracted successfully. 

Earlier interest in 7r° — * e + e~ was due to experimental 
indications 0, [nj that the decay rate could be sub- 
stantially higher than predicted, indicating possible new 
physics. Later experiments [eH fl3L Il4| obtained results 
more consistent with the standard model predictions, and 
the most recent result from KTeV-E799 Q provided a 
precise measurement of the branching ratio falling en- 
tirely within the standard model prediction. 



II. THE KTEV-E799 EXPERIMENT 

The KTeV facility (Figure at Fermilab was a 
general-purpose neutral kaon beam and decay spectrom- 
eter. It was operated for two experiments: E832, which 
used an active regenerator to produce a K$ flux for mea- 
suring SR(e'/e), and E799-II, which had a higher intensity 
Kl flux and performed searches and measurements for a 
variety of rare Kl decays. The analysis described here 
uses E799-II data. 



A. The kaon beam 

At KTeV, 800 GeV protons hit a BeO target and pro- 
duced two nearly parallel neutral beams that were de- 
fined by sweeper magnets and collimators. A vacuum 
decay volume was located from 94 to 158 meters down- 
stream of the target. In this region, the beams consisted 
of K\ and neutrons, with small numbers of shorter-lived 
neutral hyperons and Kg remaining. The K ° energies in 
this region ranged from 20-200 GeV. The decay region 
ended at a Mylar-Kevlar vacuum window which was fol- 
lowed by a charged particle spectrometer. 



B. The charged spectrometer 

Charged tracks were detected by four drift chambers 
separated by 6 m, 9 m, and 6 m. A momentum analysis 
dipole magnet sat between the second and third cham- 
bers. The field integral from the magnet was 205 MeV/c 
in the 1997 run period and 150 MeV/c in 1999. The 
momentum resolution of the spectrometer in the range 
of interest was 0.9%. A set of transition radiation de- 
tectors (TRDs) was placed after the the last drift cham- 
ber. This detector provided particle identification used 
to distinguish electrons from pions but was not needed 
in this analysis because there were no significant non- 
electron backgrounds. Following the TRDs there was a 
segmented array of scintillator planes for fast triggering 
on events with charged particles. 



C. Photon and muon detection 

The final detector for electromagnetic particles was a 
calorimeter consisting of 3100 pure Csl crystals. The 
crystal blocks were arranged in a 1.8 x 1.8 m 2 square array 
with two 15x15 cm 2 holes near the middle for the neutral 
beams to pass through. The crystals were 27 radiation 
lengths deep, which contained nearly all electromagnetic 
showers. The energy resolution for electromagnetic parti- 
cles was (0.45©2.0/a/E)%, where E is the energy in GeV 
and the addition is in quadrature. The perimeters of the 
vacuum decay region, spectrometer, and calorimeter were 
instrumented with a total of nine lead-scintillator photon 
veto counters to reject particles escaping the detector at 
high angles. Two vetoes were also used around the edges 
of the two beam holes in the calorimeter. Downstream 
of the calorimeter, a 15 cm lead wall showered remaining 
hadrons and the showers were detected by a scintillator 
plane in order to reject events with hadrons in the final 
state. Behind an additional 4 m of steel was a muon veto 
system, which was used to detect decays with muons in 
the final state. For a more complete discussion of the 
KTeV detector see Ref. [l5|. 



III. DATA COLLECTION AND ANALYSIS 

The 7T° — ► e + e~ branching ratio was found by normal- 
izing signal candidates to 7r° — > e + e~7 (Dalitz decays) 
with m e + e - > 65 MeV/c 2 . Both samples were from 
Kl 37T° decays where the other two 7r°'s in the event 
decayed to 77. This normalization mode was selected be- 
cause its final state particles and kinematics are similar 
to the signal, allowing many detector response systematic 
effects to cancel. 



A. Trigger requirements 

The trigger for both signal and normalization required 
activity in the chambers consistent with two tracks, plus 
total energy in the calorimeter above 25 GeV and at least 
four separate energy clusters in the calorimeter where at 
least one crystal in the cluster had more than 1 GeV of 
energy. The trigger also required no significant energy 
in either the photon veto counters or the hadron anti. 
Signal and normalization candidates were collected, re- 
constructed, and analyzed in parallel. 



B. Radiative corrections 

The presence of internal bremsstrahlung off the elec- 
trons in 7T° — ► e + e~ complicates the analysis, because the 
final state contains the same particles as the tree-level 
Dalitz decay 7r° — > e + e~j (though the two decays gen- 
erally populate different regions of phase space). The 
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FIG. 1: The KTeV spectrometer in E799-II configuration. 



signal must therefore be denned as a region where ra- 
diation is soft and where there is little contribution 
from the Dalitz decay. Following the conventions of 
Ref. 0, H, we defined the signal by requiring Dalitz 
Xr> = (m e + e -/m w o) 2 > 0.95, considering the rest of 
the spectrum as background. This definition left very 
little intrinsic background from the Dalitz decays while 
including 89% of the 7r° — * e + e~ bremsstrahlung spec- 
trum. Also, in this region the quantum mechanical in- 
terference between the two modes is ncglig ible 0. In 
addition to the inner bremsstrahlung diagram, a virtual 
photon correction suppresses the total ir° — > e + e~ decay 
rate by 3.4%. Both effects must be accounted for in com- 
paring the measured decay rate with theoretical models 
that neglect radiation. 

The experimentally measured quantity was the ratio: 

r(7r°^e+e-, x D >0.95) 
r(7r° e+e- 7 , x D > 0.232) ^ ' 

where the 7r° — ► e + e~ 7 rate is inclusive of ir° — ► e + e~77 
as calculated from lowest-order radiative corrections [Ta| . 



C. Reconstruction and event selection 

For both modes the full Kl — > Sir decay chain was re- 
constructed. Signal events had 6 electromagnetic clusters 
and 2 oppositely charged tracks, while the normalization 
had 7 clusters and 2 tracks. The tracks in both modes 
also had to be electron candidates, defined to be the case 



when a track of momentum p pointed to a calorimeter 
cluster of energy E and \E/p — 1| < 0.08. The total en- 
ergy in the calorimeter was required to be above 35 GeV 
and each cluster energy above 1.75 GeV. 

D. Photon and vertex reconstruction 

Clusters with no tracks pointing to them were assumed 
to be photons coming from ir° decays. For ir° — ► e + e~ 
candidates, the four photons could be assigned in three 
possible pairing combinations, while for tt° — > e + e~7 
candidates there were 15 pairing combinations for the 
five photons. The best pairing was found using the fol- 
lowing procedure: For each pair of photons the distance 
d from the calorimeter to the decay vertex was calcu- 
lated assuming the pair originated from a ir° — > 77 decay: 
d = T\2\[E\Ei 1 'm^o , where r\% was the distance between 
the two photon clusters and E\ and E% were the cluster- 
energies. The ^-position of the decay vertex was then 
z = zesi — d. A pairing \ 2 was calculated for the hypoth- 
esis that the two decay positions {z\ and Z2) coincided 
with each other and with the decay vertex obtained from 
the electron tracks (z ee ): 

x 2 = ( Z l - zf + (jg - Z? + (Zee - zf 
<7 2 {zi) CT 2 (z 2 ) <J 2 {z e e) 

For each pairing case, the mean decay position z was 
found by minimizing the x 2 ■ The pairing with the small- 
est minimum \ 2 was selected and the obtained decay 
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FIG. 2: Dalitz xz> = (m e + e - /m w o) for fully-reconstructed 
Dalitz decay candidates after normalization mode analysis 
cuts. Monte Carlo prediction using the world average form 
factor slope parameter is overlaid. 

vertex z-position, z, was then used to reconstruct parti- 
cle trajectories. This decay vertex calculation combined 
information from the calorimeter and drift chambers to 
optimize the overall resolution on the vertex position. 
The vertex was required to be 96 < z < 158 m down- 
stream of the target, removng events near the ends of 
the decay region. 

E. Final sample selection 

For 7T° — * e + e~ candidates, the reconstructed kaon 
mass was required to be between 490-510 MeV/c 2 . 
For normalization ir° — > e + e~7 candidates, where back- 
grounds were low and event reconstruction was poorer 
due to the additional pairing ambiguity, the allowed in- 
terval was 475-525 MeV/c 2 . The total reconstructed mo- 
mentum transverse to the incident kaon direction, defined 
as the line between the center of the target and the de- 
cay vertex, was required to be p\ < 10~ 3 GeV 2 /c 2 . For 
the normalization sample the reconstructed Dalitz de- 
cay mass TO e + e - 7 was required to be in the interval 100- 
200 MeV/c 2 , with an additional requirement that the 
reconstructed electron pair mass m e + e - be greater than 
70 MeV/c 2 . This last requirement removed resolution 
effects near the 65 MeV/c 2 cutoff. 

A detailed description of the detector and beamlinc 
was implemented in a Monte Carlo (MC) simulation, 
which was used to study detector geometry, acceptance, 
and backgrounds. The decay simulation included 0{a) 
radiative corrections to ir° — > e + e~ based on the work of 
Bergstrom 4], while for 7r° — > e + e~7, the world average 
form factor slope 0] and radiative corrections to order 
0(a 2 ) |l6j were used. Figure [21 shows the distribution 
of xb for the normalization Dalitz-decay sample in both 
data and MC. 

Beyond the basic reconstruction requirements above, 
additional cuts were needed to remove backgrounds to 
7T° — > e + e~. The full reconstruction of the Kl decay 



chain removed all significant backgrounds except those 
originating from Kl — > 3ir° decays. One category of 
backgrounds was Kl — > 3ir° decays with four electrons in 
the final state, where two were lost and the remaining two 
mimicked the ir° —>■ e + e~ decay. Low-energy electrons 
could be swept out of the fiducial region by the analy- 
sis magnet, never making a complete track. One major 
source of this background was Kl — > 3tt° — > (e + e~j) + 
(e + e~7) + (77), where the photons from the Dalitz de- 
cays were accidentally reconstructed as a ir° — > 77 decay. 
Another source was the rarer decay 7r° — > e + e~e + e~. Fi- 
nally, photons from 7r° decays could convert to e + e~ pairs 
in the vacuum window just upstream of the chambers. 
Events with two of these conversions, or one in combina- 
tion with a Dalitz decay, also contributed to four-track 
background. 

Backgrounds in which the two electrons came from dif- 
ferent 7r°s were reduced by a requirement on the pairing 
X 2 defined above: a cut of x 2 < 20 was used in both the 
signal and the normalization mode. To reduce the four- 
electron backgrounds further, a cut on evidence for extra 
in-time activity in the second drift chamber was made. 
Removing events with in-time activity in the second drift 
chamber more than 0.5 cm away from any reconstructed 
track reduced the four-track backgrounds to 0.7% of the 
expected signal. The effect of this cut on the signal (and 
any backgrounds without extra charged particles) was an 
overall reduction of 7.7% in acceptance. Both of these 
background cuts were also used in the Dalitz normaliza- 
tion sample in order to cancel systematic effects associ- 
ated with modeling the cut efficiency. 



F. Background and systematic error estimation 

After all cuts were applied, the largest remaining back- 
ground came from high m e + e - Dalitz decays where the 
Dalitz photon was lost and the e + e~ mass was recon- 
structed 0-0.5 MeV/c 2 high. 

A plot of m e + e - after all cuts, Fig. shows the 
signal peak at the 7r° mass and a background distri- 
bution that extends under the peak. The background 
MC normalized by the measured number of sideband 
Dalitz decays is plotted as well. The signal region was 
131.6 < m e+e - < 138.4 MeV/c 2 , in which 794 events 
were found. The MC predicted a 2.94% detector accep- 
tance for the signal in the 1997 run period and 3.14% in 
1999. In the normalization sample, 1 874 637 candidates 
were found with 0.1% background. The acceptance for 
the normalization was 1.21% in 1997 and 1.38% in 1999. 
The background in the signal region was estimated using 
a MC simulation of each of the considered backgrounds. 
Of these background events, 79% were high e + e~-mass 
Dalitz decays and the remainder were four-electron final 
states. 

The important systematic error sources that were iden- 
tified are listed in Table |U External systematic errors 
are separated so the result may be corrected in the fu- 
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FIG. 4: Ratio of data to MC distribution of m e + e - in the 
sideband region below the signal peak, where the MC was nor- 
malized to the number of fully-reconstructed Dalitz decays. 
The dotted line indicates the ratio over the entire region. 



FIG. 3: Positron-electron invariant mass for 7r° — > e + e~ sig- 
nal candidates passing all other cuts. The points with error 
bars are data; the solid histogram is background MC. 



Branching ratio uncertainties 




Statistical uncertainty 


3.8% 


7T° — > e + e _ 7 branching ratio 


2.7% 


7T° slope parameter 


1.3% 


Total external systematic uncertainty 


3.0% 


Background normalization 


1.1% 


m e + e - resolution 


0.7% 


Photon pairing \ 2 modeling 


0.5% 


Kaon momentum spectrum 


0.4% 


m e + e - cutoff in normalization 


0.3% 


Background MC statistics 


0.4% 


Signal/normalization MC statistics 


0.3% 


Total internal systematic uncertainty 


1.6% 


Total systematic uncertainty 


3.4% 


Total uncertainty on S(ty° — > e + e~) 


5.1% 



TABLE I: List of uncertainties in the ir — » e + e branching 
ratio. 



ture if the branching ratio of the Dalitz decay and the 
fraction of the decay in the high-xu region of phase 
space are measured more precisely. The Dalitz branch- 
ing ratio used was B{n° -> e+ e -j) = (1.198 ± 0.032)% 
where the relative error, 2.7%, transfers directly into 
the 7T° — > e + e~ branching ratio. The MC based on 
Ref. ^(| was used to determine the fraction of Dalitz 
events that had m e + e - > 65 MeV/c 2 , and this number 
depended on the tt° form factor used. The result was 
r(m e+e - > 65 MeV/c 2 )/T(all Dalitz) = 0.0319 when 
using the 2004 PDG[17j average for the tt° form factor 
slope. The slope value is dominated by a measurement 
in a region of spacelike momentum transfer [Tsj where an 
extrapolation using vector meson dominance was done. 
Our observed m e + e - distribution disagreed with MC at 



the 1.8a level and indicated a value that would change 
the fraction of events in the m e + e - > 65 MeV/c 2 tail by 
1.3%. This disagreement is quoted here as a systematic 
error. The detector acceptance depended negligibly on 
the form factor. 

The remaining systematic errors were internal to the 
experiment. The combination of charged and neutral in- 
formation in calculating the decay vertex caused a small 
shift in the m e + e - distribution, with the data moving 
by 0.2 MeV/c 2 more than the MC. The signal region in 
data was shifted accordingly to compensate, and an un- 
certainty in the signal acceptance and the background 
estimate was a consequence. The shift changed the ac- 
ceptance by 0.4% and the background estimate by 10.9%. 
The two errors combined into a 0.7% bias on the branch- 
ing ratio, which was taken as a systematic error. 

A systematic error was associated with the choice of 
normalization for the background. Normalizing the pre- 
diction to the number of fully-reconstructed Dalitz de- 
cays resulted in an estimate of 44.4 ± 2.7 background 
events in the signal region, where the error is from MC 
statistics only. However, the data indicated a clear ex- 
cess of events in the sideband region, 110 < m e + e - < 
130 MeV/c 2 , over this Dalitz-normalized MC. The over- 
all level of background had to be scaled up by a factor 
of 1.19 ± 0.04 (1.24 in the 1997 data; 1.15 in the 1999 
data) to match the data. The relative excess showed lit- 
tle m e + e - dependence (see Fig. 0J. There was no ex- 
cess of events in the sideband above the signal peak, 
which might have indicated an unsimulated flat contin- 
uum background. The source of the low-m e + e - sideband 
excess was not fully understood, but was likely related to 
modeling of the sensitivity of the veto system and Csl to 
the soft photon from high-xu Dalitz decays. The entire 
shift was taken as a conservative systematic error. This 
contributed a 1.1% systematic uncertainty to the branch- 
ing ratio. The final background estimate was 52. 7± 11.2. 

The high tail of the pairing \ 2 distribution was not 
simulated perfectly in the normalization and was a source 
of systematic uncertainty. Removing the x 2 cu t m t ne 
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normalization analysis changed the measured number of 
decaying kaons by 0.5%. This was not expected to cancel 
in the ratio, as the pairing \ 2 distributions were different 
between signal and normalization due to the presence of 
an additional photon in the normalization sample. The 
entire sensitivity of the normalization level to the cut was 
taken as a systematic error. 

The simulated kaon momentum distribution deviated 
from the data, as evidenced by a slope in the ratio of 
the reconstructed momentum distributions in data and 
MC. Each MC event was reweighted to account for the 
slope in both signal and normalization. This modification 
changed the ratio of signal to normalization acceptances 
by 0.4%, which was taken as a systematic error on the 
branching ratio. In the normalization, the cut on m e + e - 
caused a small bias in the branching ratio due to mod- 
eling of the acceptance near the m e + e - = 70 MeV/c 2 
boundary. Tightening the cut by 5 MeV/c 2 produced a 
0.4% difference in the branching ratio. 



of xd yields 

B(tt° -> e+e~, x D > 0.95) = (6.44 ± 0.25 ± 0.22) x 10~ 8 

(4) 

In both cases the first error is from data statistics alone 
and the second is the total systematic error. 

Comparison with theoretical predictions and the uni- 
tary bound can be done only if we remove the effects of 
final state radiation. This was done by extrapolating the 
full radiative tail beyond xd — 0.95 and scaling the re- 
sult back up by the overall radiative correction of 3.4% 
to find the lowest-order rate for tt° — ► e + e _ . We found 
gno-rad^o _^ e + e -) = (7.48 ± 0.29 ± 0.25) x 10~ 8 , more 
than seven standard deviations higher than the unitary 
bound. The result falls between VMD 6] and xPT pre- 
dictions Q, with a significance on the difference of 2.3 
and 1.5 standard deviations respectively. 



IV. RESULTS AND CONCLUSION 
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The final branching ratio was calculated from 794 can- 
didate signal events with an estimated background of 
52. 7± 11.2, and 1 874 637 normalization events with neg- 
ligible background. We found 

ry^« + «- ^> 0-95) = (1 . 685± o.064±0.027)xlO- 
r(7r° -> e+e-j,x D > 0.232) v ' 

(3) 

where Xd — 0.232 corresponds to m e + e - = 65 MeV/c 2 . 
Extrapolating the Dalitz branching ratio to the full range 
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